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(57) The demodulator of the invention includes a 
DSP. The DSP determines the distortion amount of the 
transmission channel based on known signals included 
in a radio receiving signal, and estimates the frequency 
of the direct wave based on the distortion amount. Fur- 
thermore, the DSP rotates the phase of the distortion 
amount based on the phase amount corresponding to 
the estimated frequency of the direct wave. Next, the 
DSP estimates the center frequency of the Doppler 
spread based on the distortion amount after the phase 
is rotated. Furthermore, the DSP synthesizes the phase 
amount corresponding to the estimated center frequen- 
cy of the Doppler spread and the phase amount corre- 
sponding to the frequency of the direct wave to estimate 
the last frequency offset. Thereafter, the DSP eliminates 
frequency offset of the radio receiving signal based on 
the estimated frequency offset. Thereby, a sufficient fre- 
quency offset compensating range can be secured, and 
excellent BER characteristics can be realized. 



FIG.1 



z 
o 

to I- 

625 

2 O 

co <r 

I 5 



T 



z>tr 

8° 



z cr 
2o 
"z 

p 

Z UJ 

to 
z 




Printed by Jouve, 75001 PARIS (FR) 



BNSDOCID: <EP 1180853A1J_> 



10 



15 



35 



40 



45 



50 



55 



EP 1 180 853 A1 

Description 

BACKGROUND OF THE INVENTION 
1 .Field of the Invention 



[0001] The present invention relates to an automatic frequency control method and device, which are applied to 
digital demodulation processing in satellite communication, mobile satellite communication, mobile terrestrial commu- 
nication, and a demodulator to which said device is applied. 

2. Description of the Related Art 



[0002] Recently, in satellite communication, mobile satellite communication, and mobile terrestrial communication, 
digital modulation and demodulation have been actively researched. In particular, in the environment of mobile com- 
munications, a radio signal is greatly influenced by fading. Therefore, various demodulation systems which stably 
operate even in such a fading environment have been considered. Among these systems, a system, which is con- 
structed so that absolute coherent detection can be executed even in a fading environment by estimating and com- 
pensating fading distortion by using known signals, has been recognized. In the case where fading distortion is esti- 
mated and compensated after quasi-coherent detection or the like is carried out by this system, to estimate and com- 
20 pensate fading distortion with high accuracy, it is necessary that the frequency offset between the carrier wave fre- 
quency of a radio transmission signal and the oscillation frequency of a reference signal for quasi-coherent detection 
is small. 

[0003] However, in the case where the frequency stability and accuracy of the oscillation circuit of a transmitter- 
receiver are insufficient, unless the frequency of a radio receiving signal is automatically controlled by eliminating this 
25 frequency offset by some processing, fading distortion cannot be estimated and compensated with high accuracy. 

[0004] In mobile communication, transmission and receiving are carried out between a fixed station and a mobile 
station or between mobile stations. Therefore, when two stations relatively move, the frequency of a radio receiving 
signal deviates due to Doppler variation. Therefore, even if the stability and accuracy of the oscillation circuit of the 
transmitter-receiver are good, frequency offset occurs between the frequency of the radio receiving signal and the 
so oscillation frequency of the reference signal. 

[0005] A technique for compensating frequency offset is disclosed in, for example, Japanese Patent Laid-Open No 
93302/1 997 titled "DIGITAL MOBILE RADIO COMMUNICATION METHOD" . By the conventional technique disclosed 
in this document, frequency offset is eliminated by using phase fluctuation information of known signals (pilot signals). 
[0006] In this conventional technique, a radio transmission signal into which two-symbol known signals have been 
inserted every inserting period N F from the transmission side is transmitted. On the other hand, at the receiving side, 
the phase change amount between the sequential two-symbol known signals is calculated, and in accordance with 
the calculated phase change amount, the phase of the radio receiving signal is rotated. Frequency offset is thus elim- 
inated from the radio receiving signal. 

[0007] In electric wave transmission channels, there is a Rician fading transmission channel in which direct waves 
and multi-path waves are mixed. In this case, the direct wave is Doppler-shifted. Therefore, the frequency f D of the 
direct wave further deviates by Doppler-shift amount f DP caused by the Doppler-shifting from the offset f OFST caused 
by stability of the oscillation circuit as shown in Fig. 24A. 

[0008] On the other hand, in the abovementioned prior art, the phase change amount between the sequential two- 
symbol known signals is supposed as frequency offset. The frequency offset calculated in this case is equivalent to 
the difference between the oscillation frequency f G at the receiver side and the frequency f D of the direct wave. That 
is, in the prior art, as shown in Fig. 24B, frequency control is carried out so that the frequency f D of the direct wave 
almost coincides with the oscillation frequency f c . In this case, the center frequency f M of the Doppler spread deviates 
from the oscillation frequency f c by the Doppler-shift amount f DP . Therefore, the Doppler spread apparently further 
broadens, so that the frequency corresponding to the end of the Doppler spread greatly deviates from the oscillation 
frequency f 0 . Accordingly, frequency offset cannot be satisfactorily compensated. Therefore, bit error rate character- 
istics (hereinafter, referred to as BER characteristics) after the radio receiving signal is demodulated may deteriorate. 

SUMMARY OF THE INVENTION 



[0009] The object of the invention is therefore to solve the abovementioned technical problems and provide an au- 
tomatic frequency control method and device by which satisfactory BER characteristics can be realized without being 
influenced by a Doppler-shifted direct wave in a Rician fading environment where the direct wave has Doppler-shifted. 
[001 0] Another object of the invention is to provide a demodulator by which demodulation accuracy can be improved 
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by applying the abovementioned automatic frequency control device to this demodulator. 

[0011] In order to achieve the abovementioned objects, an automatic frequency control method of the invention, 
which controls the frequency of a radio receiving signal by compensating frequency offset of the radio receiving signal 
that periodically includes a plurality of adjacent known signals, comprising, based on the distortion amount of each 
5 known signal included in the radio receiving signal, the frequency of the direct wave of the radio receiving signal and 
the center frequency of the Doppler spread of the radio receiving signal are estimated, and based on both of the 
frequencies, frequency offset of the radio receiving signal is compensated. 

[0012] In this construction, by compensating the frequency offset based on the frequency of the direct wave, the 
compensating range of frequency offset can be sufficiently secured, and by compensating the frequency offset based 
10 on the center frequency of the Doppler spread, satisfactory BER characteristics can be secured. That is, compatibility 
with securing of a sufficient frequency offset compensating range and securing of satisfactory BER characteristics can 
be realized. 

[0013] In addition, the abovementioned construction is applied to a demodulator, whereby demodulation processing 
can be executed for a radio receiving signal in which frequency offset has been satisfactorily eliminated. Therefore, 
is the demodulation quality can be improved. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] 

20 

Fig. 1 is a conception diagram showing the entire structure of the radio communication system to which an auto- 
matic frequency control method relating to Embodiment 1 of the invention is applied; 
Fig. 2 is a drawing showing the format of a burst signal; 
Fig. 3 is a flowchart for explanation of demodulation processing; 
25 Fig. 4 is a drawing for explanation of estimated timing of a frequency offset; 

Fig. 5 is a flowchart for explanation of a frequency offset estimation processing; 
Fig. 6 is a graph showing BER characteristics; 

Fig. 7 is a flowchart for explanation of the first frequency offset estimation processing; 
Fig. 8 is a flowchart for explanation of the second frequency offset estimation processing; 
30 Fig. 9 is a conception diagram for explanation of the second frequency offset estimation processing; 

Fig. 1 0 is a graph for explanation of BER based on Embodiment 1 ; 

Fig. 1 1 is a flowchart for explanation of the first frequency offset estimation processing relating to Embodiment 2 
of the invention; 

Fig. 1 2 is a flowchart for explanation of the second frequency offset estimation processing relating to Embodiment 
35 3 of the invention; 

Fig. 1 3 is a flowchart for explanation of the second frequency offset estimation processing relating to Embodiment 
4 of the invention; 

Fig. 14 is a flowchart for explanation of the first frequency offset estimation processing relating to Embodiment 5 
of the invention; 

40 Fig. 15 is a flowchart for explanation of the first frequency offset estimation processing relating to Embodiment 6 

of the invention; 

Fig. 16 is a flowchart for explanation of frequency offset estimation processing relating to Embodiment 7 of the 
invention; 

Fig. 17 is a flowchart for explanation of demodulation processing relating to Embodiment 8 of the invention; 
45 Fig. 1 8 is a flowchart for explanation of demodulation processing relating to Embodiment 9 of the invention; 

Fig. 1 9 is a flowchart for explanation of demodulation processing relating to Embodiment 1 0 of the invention; 
Fig. 20 is a flowchart for explanation of the first AFC processing; 
Fig. 21 is a flowchart for explanation of the second AFC processing; 

Fig. 22 is a flowchart for explanation of demodulation processing relating to Embodiment 11 of the invention; 
so Fig. 23 is a flowchart for explanation of demodulation processing relating to Embodiment 12 of the invention; and 

Fig. 24 is a conception diagram for explanation of the relationship of the frequency offset estimating range, esti- 
mation accuracy, and possible frequency offset. 

DETAILED DESCRIPTION OF THE PREFFERED EMBODIMENTS 

55 

[0015] Hereinafter, Embodiments of the invention are explained in detail with reference to the attached drawings. 
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Embodiment 1 

[0016] Fig. 1 is a block diagram showing the construction of a radio communication system to which an automatic 
frequency control method relating to Embodiment 1 of the invention is applied. This radio communication system com- 
prises transmitter 1 and receiver 1 0, and has a function for estimating and eliminating frequency offset of a burst signal 
radio-transmitted from the transmitter 1 by receiver 10. 

[0017] Concretely, for this radio communication system, a satellite communication system, mobile satellite commu- 
nication system, and mobile terrestrial communication system can be applied. In the satellite communication system, 
transmitter 1 and receiver 1 0 are applied to earth stations set at different positions on the earth. In the mobile satellite 
communication system, the transmitter 1 is applied to either an earth station or a mobile station that is set on the earth, 
and the receiver 10 is applied to an opposite station of the transmitter 1 of the two stations. In the mobile terrestrial 
communication system, the transmitter 1 is applied to either a base station or a mobile station, and the receiver 10 is 
applied to an opposite station of the transmitter 1 of the two stations. 

[0018] As a mobile station in the mobile satellite communication system, a satellite mobile phone of a single mode 
or a dual mode can be applied. In addition, a cellular phone can be applied to the mobile station in the mobile terrestrial 
communication system. 

[0019] This communication system uses the TDMA (Time Division Multiple Access) as a communication method. 
Therefore, the transmitter 1 radio-transmits a burst signal synchronized with a predetermined time slot to the receiver 
10. The receiver 10 demodulates the receiving burst signal, eliminates frequency offset, and then restores the original 
signal. 

[0020] The transmitter 1 comprises known signal insertion circuit 2, modulation circuit 3, and transmission circuit 4. 
The known signal insertion circuit 2 periodically inserts known signals into information signals to be transmitted, and 
generates burst signals before modulation. Concretely as shown in Fig. 2A, the known signal insertion circuit 2 inserts 
Np-symbol sequential known signal rows (hereinafter, referred to as known signal blocks) into (N F -N p )-symbol infor- 
mation at every N F -symbol periods. N P is an integer of 2 or more (N P >2). The timings at which the known signals are 
inserted are shown by (kNp+i)T s . Herein, k indicates the insertion order of the known signals, i is a value between zero 
and (Np-1) (0<i<N P -1). T s is a symbol period. 

[0021] The known signal insertion circuit 2 provides a burst signal before being modulated to the modulation circuit 
3. The modulation circuit 3 modulates this non-modulated burst signal, and outputs it as a burst signal after being 
modulated. The modulation circuit 3 synchronizes the modulated burst signal with a predetermined time slot and pro- 
vides the signal to the transmission circuit 4. Concretely, as shown in Fig. 2B the modulation circuit 3 synchronizes 
burst signals B. B+1. B+2, B+3... with predetermined time slots S1 5 S2, S3, S4... and provides the signals to the 
transmission circuit 4. The transmission circuit 4 superposes the burst signals on electric waves and transmits them 
to the receiver 1 0. 

35 [0022] The receiver 10 includes receiving circuit 11 and demodulator 12. The receiving circuit 11 is comprised of an 
amplifier, a frequency converter, and the like, and converts the frequency of the receiving burst signal which is a radio 
receiving signal into an intermediate frequency and outputs it as a receiving I F signal. The demodulator 1 2 demodulates 
the receiving IF signal outputted from the receiving circuit 11 by quasi-coherent detection, and restores the original 
information signal. More concretely, the demodulator 12 comprises frequency converter circuit 21 , two A/D (analog/ 
digital) converter circuits 25a and 25b, and DSP (Digital Signal Processor) 26 as an automatic frequency control device 
or digital signal processing device, and restores the original information signal by digitally executing demodulation 
processing by the DSP 26. 

[0023] The frequency converter circuit 21 converts a receiving IF signal periodically including a plurality of known 
signals into an analog baseband signal, and the circuit includes one oscillation circuit 22, two multiplication circuits 
23a and 23b, and n/2 phase shift circuit 24. The oscillation circuit 22 generates local oscillation signals with a prede- 
termined oscillation frequency. The local oscillation signals generated by the oscillation circuit 22 are inputted into two 
multiplication circuits 23a and 23b. In this case, the local oscillation signal is inputted into the multiplication circuit 23a 
for I channel via the n/2 phase shift circuit 24. The n/2 phase shift circuit 24 shifts the phases of the local oscillation 
signals by n/2. Therefore, the local oscillation signals whose phases are shifted by n/2 from each other are inputted 
50 mio the multiplication circuits 23a and 23b. 

[0024] The multiplication circuits 23a and 23b each mix the receiving IF signals and the local oscillation signals. As 
a result, I channel and Q channel analog baseband signals are generated. These generated analog baseband signals 
are supplied to the A/D converter circuits 25a and 25b respectively. The A/D converter circuits 25a and 25b each 
convert the analog baseband signals into I channel and Q channel digital baseband signals. The digital baseband 
55 signals are supplied to DSP 26 respectively. 

[0025] The DSP 26 receives the inputs of the digital baseband signals generated by the A/D converter circuits 25a 
and 25b, and based on the distortion amounts of the known signals included in the inputted digital baseband signals 
estimates the frequency of the direct wave and the center frequency of the Doppler spread, and compensates frequency 
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offset based on these estimated frequencies. Thus, the DSP 26 automatically controls the frequencies of the digital 
baseband signals. Furthermore, the DSP 26 eliminates fading distortion from the digital baseband signals from which 
frequency offset has been compensated, and then demodulates the digital baseband signals, whereby the original 
information signals are restored. 

5 [0026] More concretely, the DSP 26 has storage 26a comprised of a ROM or the like. In the storage 26a, a demod- 
ulation processing program which is a computer program is stored. By executing the demodulation processing program 
stored in the storage 26a, the DSP 26 executes serial demodulation processing including estimation of the abovemen- 
tioned two frequencies, compensation for frequency offset, compensation for fading distortion and demodulation. 
[0027] Fig. 3 is a flowchart for explanation of the demodulation processing to be executed by the DSP 26. The DSP 

10 26 applies filtering such as waveform shaping processing to the digital baseband signals supplied from the A/D con- 
verter circuits 25a and 25b (step S1). Thereby, the DSP 26 eliminates the high frequency components such as noise 
components higher than a predetermined cutoff frequency from the digital baseband signals. 

[0028] Thereafter, the DSP 26 executes the Nyquist point detection processing (step S2). Concretely, the DSP 26 
detects the Nyquist point of the digital baseband signal to obtain digital baseband signal r(kN F +i) corresponding to the 
15 Nyquist point. The digital baseband signal r(kNp+i) corresponding to the Nyquist point is expressed as Formula (1) 
below. In Formula (1) below, c(kNp+i) shows a distortion amount caused by fading. In addition, A shows the amplitudes 
of the signals, and b(kNp+i) shows symbol values. Furthermore, n(kNfH-i) shows gauss noises. 

2Q r(kN F + i) = Ac(kN F + i)b(kN F + i) + n(kN F + i) (1 ) 

[0029] Thereafter, the DSP 26 executes automatic frequency control processing (Step S3 through S5). In the auto- 
matic frequency control processing, frequency offset based on the oscillation frequency f 0 of the oscillation circuit 22 
as a reference is estimated, and the estimated frequency offset is eliminated from the digital baseband signal r(kNp-i-i), 

25 whereby the frequency of the digital baseband signal r(kNp+i) is automatically controlled. In other words, the DSP 26 
estimates the frequency of the direct wave and the center frequency of the Doppler spread based on the distortion 
amounts of the plurality of known signals periodically included in the digital baseband signal corresponding to the radio 
receiving signal, compensates the frequency offset of the digital baseband signal r(kNp+i) based on both the above- 
mentioned frequency, and thereby automatically controls the frequencies of the digital baseband signal r(kNpH). 

30 [0030] Further in other words, the DSP 26 estimates frequency offset of the digital baseband signal r(kNp+i) from 
the frequency of the direct wave and the center frequency of the Doppler spread that are estimated based on the 
distortion amounts of the plurality of known signals periodically included in the digital baseband signal r(kN F +i) corre- 
sponding to a radio receiving signal, eliminates the estimated f req uency offset from the digital baseband signal rfkNp+i) , 
whereby the frequency of the digital baseband signal r(kNp+i) is automatically controlled. 

35 [0031] More concretely explaining the automatic frequency control processing, the DSP 26 executes frequency offset 
estimation processing first (step S3). More specifically, the DSP 26 estimates a phase difference e s (ml_N F ) correspond- 
ing to the phase rotation amount of one symbol as frequency offset based on the digital baseband signal rfkNp+i). Still 
more specifically, the DSP 26 determines a phase difference 0 s (mLN F ) as frequency offset at each estimating time mT 
(m: natural number) that comes after each estimating period T (=LN F T S ) (see Fig. 4) based on the distortion amounts 

40 of the plurality of known signals included in the digital baseband signal rfkNjr+i). 

[0032] Next, the DSP 26 executes integration processing for integrating the determined phase difference G s (mLN F ) 
by one-symbol period T s (step S4). Concretely, the DSP 26 cyclically sums phase differences 6 s (mLN F ) at each symbol 
period T s as shown in Formula (2) below. Thereby, the DSP 26 obtains an accumulated phase difference 6((m-1) 
LN F +1 Np+i). In Formula (2) below, 1 indicates identification numbers corresponding one to one known signal blocks, 

45 and is a value of 0 or more and (L-1) or less (0<1<L-1). L shows the number of known signal blocks included in the 
estimating period T. 

8((m - 1)LN F + 1N F + i) = 8((m - 1)LN F + t N F + i - 1) + 6 s (mLN F ) (2) 

50 

[0033] Thereafter, the DSP 26 executes frequency offset elimination processing for eliminating frequency offset from 
the digital baseband signal r(kN F +i) based on this accumulated phase difference e((m-1)LN F +IN F +i) (step S5). Con- 
cretely, the DSP 26 rotates in reverse the phase of the digital baseband signal r(kNp+i) by an amount in accordance 
with the accumulated phase difference 6((m-1 )LN F +1 N F +i). Thereby, frequency offset can be eliminated from the digital 
55 baseband signal r(kNp+i). That is, as shown in Formula (3) below, the DSP 26 obtains a digital baseband signal r R 
(kNp+i) in which frequency offset has been eliminated. Thus, the DSP 26 automatically controls the frequency of the 
digital baseband signal. 
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r^(kN F + i) = r K ((m - lJLN r + IN, + i) 

- r< (m - l)LN p + 1N F + i) exp[-j8( (m - 1>LN P + lN p + i) ] 

(3) 

[0034] Thereafter, the DSP 26 executes fading distortion compensation processing for estimating and eliminating 
fading distortion from the digital baseband signal r R (klVi) (step SB). Concretely, the DSP 26 extracts N P -symbol 
known signals from the digital baseband signal r R (kN F+ i). The frequency offset caused by the transmitter 1 and receiver 
1 0 has been already eliminated from these extracted known signals. 

[0035] The DSP 26 detects fading distortion based on the extracted N P -symbol known signals. The DSP 26 executes 
interpolation processmg such as gauss interpolation, Wiener interpolation, based on the detected fading distortion 
The DSP 26 thus estimates fading distortion in the information signal. Furthermore, the DSP 26 eliminatesthe estimated 

dl !! ort,on - There °y. the DSP 26 executes fading interpolation for the digital baseband signal r R (kNp+i) 
[0036] Thereafter, the DSP 26 executes data decision processing (step S7). Concretely, the DSP 26 dicides the 
original information signal based on the digital baseband signal in which fading has been compensated Thus the 
demodulation processing is completed. ' 
*» [0037] Fig. 5 is a flowchart for explanation of the frequency offset estimation processing. This frequency offset esti- 
mation processing includes the first frequency offset estimation processing and the second frequency offset estimation 
processing. That is, the DSP 26 estimates frequency offset by combining the first frequency offset estimation processing 
and the second frequency offset estimation processing. 

[0038] In the first frequency offset estimation processing, the frequency of the direct wave of the radio receiving 
signal .s estimated as the first frequency offset. In the second frequency offset estimation processing, the center fre- 
quency of the Doppler spread of the radio receiving signal is estimated as the second frequency offset The DSP 26 
estimates the last frequency offset from the first and second frequency offsets estimated by these two frequency offset 
estimation processings. y 
[0039] Thus, the reason why the two frequency offset estimation processings are combined is as follows. As explained 
in the subsection "Description of the Related Art", in the case where the frequency f D of the direct wave is regarded 
as frequency offset and compensated, as shown in Fig. 24B, the Doppler spread apparently further broadens, and the 
characteristics deteriorate. On the other hand, in the case where the center frequency f M of the Doppler broadening 
is regarded as frequency offset and compensated, as shown in Fig. 24C, the direct wave has frequency deviation 
correspondmgto the Doppler shift amount f DP , however, the Doppler spread is left at the original degree of broadening 
Therefore, in this case, the deterioration of the BER characteristics is less than that in the case where the frequency 
f D of the direct wave is regarded as frequency offset and compensated. 

[0040] Fig. 6 is a graph showing the BER characteristics. As can be clearly understood from this graph at the point 
where the compensation for a frequency close to the center frequency f M of the Doppler spread suppresses the dete- 
rioration of the BER characteristics more greatly than in the case where the frequency f D of the direct wave is com- 
pensated, it is understood that the center frequency f M is the optimum frequency offset. That is, in terms of prevention 
of deterioration in the BER characteristics, it can be said that the center frequency f M of the Doppler spread being 
regarded as the frequency offset is better than the frequency f D of the direct wave being regarded as the frequency 

[0041] On the other hand, in the case where the center frequency f M of the Doppler spread is regarded as frequency 
offset and compensated, prevention of deterioration in the BER characteristics is possible, however, there is a problem 
where the frequency offset compensating range is relatively narrow. On the other hand, in the case where the frequency 
f D of the direct wave is regarded as frequency offset and compensated, a relatively wide frequency offset compensating 
range can be secured. r a 

[0042] Concretely, on supposition that the phase rotates by A6 in a certain time At, the freauencv offset Af can h e 
expressed by hormula (4) below. 
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55 [0043] When this frequency offset Af is estimated by using the phase rotation amount between the known signals 
he offset can be expressed by Formula (5) below. In Formula (5) below, A6 P shows the phase rotation amount between 
the known signals. In addition, Fysymbol/s) shows the transmission rate of the signals 
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Af 1 !g 1 f|s Aflp (5) 
27i At 2n N F p 

[0044] The detecting range of the phase rotation amount A0 P is -7i<Ape<7c, and as a result, the estimating possible 
range of the frequency offset Af can be expressed by Formula (6) below. 

- Jk < Af < ^ (6) 
2N F " 2N F 



[0045] When the frequency f D of the direct wave is estimated, since the phase rotation amount between the adjacent 
known signals is used, the f D becomes equivalent to the value obtained by Formula (6) upon substituting 1 for N F . That 
is, the frequency offset estimating range f DEri in the first frequency offset estimation processing becomes a range 

15 between -R s /2 and R s /2. 

[0046] In addition, when the center frequency f M of the Doppler spread is estimated, since the phase rotation amount 
between the known signal blocks inserted for each N F symbol is used, the frequency offset estimating range f DET2 
becomes the same as the range expressed by the abovementioned Formula (6). That is, the frequency offset estimating 
range f DET2 in lne second frequency offset estimation processing becomes the range between -R S /2N F and R S /2N F . 

20 [0047] Thus, it is understood that the frequency offset estimating range f DET1 in the first frequency offset estimation 
processing is N F times larger than that in the second frequency offset estimation processing. That is, in the first fre- 
quency offset estimation processing, the frequency offset can be estimated in the relatively wide frequency offset 
estimating range f DET v 

[0048] From the above description, in this Embodiment 1 , both of the frequency f D of the direct wave and the center 
25 frequency f M of the Doppler spread are estimated, and both of these frequencies f D and f M are used for compensation 
for the frequency offset, whereby deterioration in the BER characteristics can be prevented while securing a sufficient 
frequency offset compensating range. 

[0049] Hereinafter, the frequency offset estimation processing is explained in detail. The DSP 26 executes distortion 
amount detection processing first (step T1). In the distortion amount detection processing, the distortion amounts of 

30 the transmission channel are detected in symbol units based on the known signals. That is, as shown by Formula (7) 
below, based on the Np-symbol known signal blocks included in the digital baseband signal rfkNp+i), the DSP 26 
detects the distortion amounts c EPj (N F +i) of the transmission channel for each symbol in the known signal blocks. In 
this case, the distortion amounts CEp^kNp+i) of the transmission channel corresponding to the amplitude and the dis- 
tortion amount of the phase of the digital baseband signal rfkNp+i). In addition, in Formula (7) below, b P shows the 

35 symbol value of the known signal. 



r(kN P + i) 

40 . n(kN F + i) < 7 > 

« Ac(kN F + x) + — L 



45 [0050] Thereafter, the DSP 26 executes the first frequency offset estimation processing based on the detected dis- 
tortion amounts c^kNp+i) of the transmission channel (step T2). Concretely, the DSP 26 estimates the frequency of 
the direct wave as the first frequency offset based on the distortion amounts c EPj (kNp+i) of the transmission channel. 
More concretely, the DSP 26 estimates the phase difference 6 EP1 (mLN F ) corresponding to the phase rotation amount 
between the adjacent known signals as the first frequency offset based on the distortion amount Cep^kNp+i) between 

so the adjacent known signals among distortion amounts c^kNp+i) of the transmission channel. 

[0051] Thereafter, the DSP 26 executes frequency offset elimination processing for eliminating the first frequency 
offset from the distortion amount of each known signal based on the estimated phase difference e EP1 (mLN F ) (step T3). 
More concretely, the DSP 26 compensates the phase by rotating the phases of the distortion amounts c EP ,(kN F 4-i) of 
the transmission channel by the phase difference 6 EP1 (mLN F ) as shown by Formula (8). Thereby, the DSP 26 obtains 

55 the distortion amounts Cep-^kNIp+i) in which the first frequency offset has been eliminated. 

c EP1 (kN F + i) = c EPI (kN F + i)exp[- j6(kN F + i)] 
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e(kN F + i) = 6(kN F + i - 1 ) + e EP1 (mLN F ) (8) 

[0052] Thereafter, the DSP 26 executes the second frequency offset estimation processing (step T4). Concretely, 
the DSP 26 estimates the second frequency offset based on the distortion amounts c EP1 (kN F +i) in which the first 
frequency offset has been eliminated. More concretely, the DSP 26 estimates the phase difference e EP2 (mLN F ) cor- 
responding to the phase rotation amount for one symbol as the second frequency offset based on the distortion amount 
of the periodically inserted known signals among the distortion amounts c EP1 (kN F +i) in which the first frequency offset 
has been eliminated. 

[0053] Thereafter, the DSP 26 executes phase synthesis processing (step T5). Concretely, the DSP 26 synthesizes 
the two phase differences e EP1 (mLN F ) and e EP2 (mLN F ) as shown by Formula (9) to estimate the phase difference G s 
(rnLNp) as the last frequency offset. 



15 e s (mLN F ) = 0 EP1 (mLN F ) + 8 EP2 (mLN F ) (9) 

[0054] Fig. 7 is a flowchart for detailed explanation of the first frequency offset estimation processing. The DSP 26 
executes phase difference vector operation processing first (step U1). Concretely, the DSP 26 determines the phase 
difference vector D EP (kN F ) based on the distortion amounts C EPj (kNpfi) of the transmission channel. More concretely, 
as shown by Formula (10) below, based on the distortion amount between the adjacent known signals in an arbitrary 
known signal block among the distortion amounts C EPi (kN F +i) of the transmission channel, the DSP 26 determines the 
phase difference vector c EPI (kN F +i). In Formula (10) below, i shows the identification numbers corresponding one to 
one known signals, and are zero or larger and (N p -2) or smaller (0<i<N P >2). In addition, shows a complex conjugate. 
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H p -2 

D EP (kN r ) - 2 o C **i + i/ kN F + 1 + 1 > C EPi * (kN p + i) (10) 

[0055] Thereafter, the DSP 26 executes phase difference vector average processing (step U2). Concretely, the DSP 
26 averages the phase difference vectors D EP (kN F ) in arbitrary known signal blocks for the estimation period T to 
determine the average phase difference vector D EPA (mLN F ). 

[0056] More concretely, the DSP 26 collects the phase difference vectors D EP (kN F ) in known signal blocks for the 
estimation period T from the estimation time (m-1 )T to the estimation time mT (see Fig. 4). If the number of known 
signal blocks in the estimation period T is L, the DSP 26 obtains L of phase difference vectors D EP ((m-1)T+1N F ). 
Herein, 1 shows the identification numbers corresponding one to one known signal blocks, and is 0 or more and (L-1) 
or less (0<1<L-1). Thereafter, as shown by Formula (11) below, the DSP 26 determines the average phase difference 
vector D EPA (mLN F ) by averaging the collected L of phase difference vectors D EP ((m-1 )LN,h-IN f ). 



L-1 

D EP*(roLN r ) = 2 D ^< m - 1 J LN f + 1N F ) (11> 

1-0 



[0057] Then, the DSP 26 executes phase difference operation processing (step U3). Concretely, as shown by For- 
mula (12) below, the DSP 26 determines the phase difference 9 EP1 (mLN F ) based on the average phase difference 
vector D EPA (mLN F ). Thus, the first frequency offset according to the frequency of the direct wave is estimated within 
the frequency offset estimating range f DET1 that is wider than the frequency offset estimating range 1 DET ?. 

„ , i K i \ «1 ,m [D EPA (mLN F )] 



[0058] Fig. 8 is a flowchart for explanation of the second frequency offset estimation processing. The DSP 26 exe- 
cutes distortion amount average processing first (step V1). Concretely, the DSP 26 applies average processing to the 
distortion amounts c EP1 (kN F +i) (0<t<N P -1) of the transmission channel in the known signal blocks in which the first 
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frequency offset has been eliminated to determine the average distortion amount CepCkNp). 

£0059] More concretely, as shown by Formula (13) below, the DSP 26 sums the distortion amounts c EP1 (kN F +i) 
corresponding to each symbol in arbitrary known signal blocks, and divides the summed value by the symbol number 
N P of the known signals in the known signal blocks. Thus, the DSP 26 obtains the average distortion amount c EP (kN F ) 
5 of the transmission channel with regard to one known signal block in which noise or the like has been eliminated. 



10 



C EP (kN F ) 



iV, 

Np i~ 0 



c E pi< kN r + i) 



(13) 



[0060] Next, the DSP 26 executes interblock discrete Fourier transform (DFT) processing based on this average 
distortion amount c EP (kN F ) (step v2). Concretely, the DSP 26 applies DFT processing to the average distortion amounts 

15 c EP ((m+1)N F )(0<1<L-1) in certain L of known signal blocks within the estimation period T (=LN F T S ), and thereby de- 
termines signal electric powers P f (n) corresponding to the plurality of frequency offset candidates nAf RES . 
[0061] More concretely, the DSP 26 determines signal powers P f (n) corresponding to the plurality of frequency offset 
candidates nAf RES that are set for each predetermined estimation accuracy Af RES as shown by the circles in Fig. 9A 
within the frequency offset estimating range foET2- Herein, the frequency offset estimating range f DET2 , as mentioned 

20 above, is regulated by the known signal insertion period N F , and if the estimation accuracies Af RES are used, this range 
is a range between -MAf RES and MAf RES . M is a constant, and is approximately expressed by Formula (14) below. 



25 



M ~ 



2N F Af RES 



(14) 



30 



35 



[0062] The DSP 26 rotates the phases of the average distortion amounts C EP ((m+1)N F ) in the frequency offset 
estimating range f DET2 by the phase amounts corresponding to the respective possible frequency offset candidates 
nAf RES . Thereafter, the DSP 26 vector-synthesizes the average distortion amounts c EP (m+1)N F ) whose phases have 
been rotated. Thereby, the DSP 26 obtains signal powers P f (n) corresponding to the plurality of frequency offset can- 
didates nAf RES as shown by the group of arrows in Fig. 9B. 

[0063] The abovementioned processing can be summarized as Formula (15) below. In Formula (15) below, R s shows 
the signal transmission rate, n is a value between -(M+W) and (M+W). In addition, W is a parameter showing the 
frequency bandwidth of the frequency window described later. 



40 



J c„((n + l)N„)exp(- j 

1-0 \ 



(15) 



[0064] Thereafter, the DSP 26 executes window signal power operation processing (step V3). Concretely, the DSP 
26 determines the window signal powers E f (n) that correspond one to one frequency windows having predetermined 
frequency widths. For example, the frequency width of a frequency window is 2W times of the estimation accuracy 
« Af RES . W is set in accordance with the fading condition of the transmission channel, and is set to be, for example, two 
times of the Doppler spread. As shown by Formula (1 6) below, the DSP 26 sums the signal powers P f (n) of the frequency 
offset candidates nAf RES existing in this frequency window to determine the window signal power E f (n) corresponding 
to this frequency window. In Formula (16) below, n is a value between -M and M. 



50 

w 

E f (n) = 2 P * (n + k) < 16 > 



55 [0065] Thus, by smoothing the signal powers P f (n) by using the frequency windows, among the window signal powers 
E f (n) corresponding to the plurality of frequency offset candidates nAf RES , as shown in Fig. 9C, the window signal 
power E f (n) of the frequency offset candidate nAf RES accordingto the center frequency of the Doppler spread becomes 
maximum. 
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[0066] Next, the DSP 26 executes maximum value detection processing (step V4). Concretely, as shown by Formula 
(17) below, the DSP 26 detects the maximum value among all window signal powers E f (n), and estimates the frequency 
offset candidate nAf RES corresponding to the maximum value as the second frequency offset. That is, the DSP 26 
estimates the center frequency of the Doppler spread as the second frequency offset. 

^OFST ~ n WAX^^lES 

E f (rW = Max[E £ (n)] (17) 



[0067] Next, the DSP 26 executes phase difference operation processing (step V6). Concretely, as shown by Formula 
(18) below, the DSP 26 determines the phase difference e EP2 (ml_N F ) of one symbol based on the estimated second 
frequency offset. Thus, the phase difference G EP2 (mLN F ) corresponding to the center frequency of the Doppler spread 
is is obtained. 

e EP2 (mLN F ) = f OFST x p (18) 

[0068] As described above, according to this Embodiment 1 , by estimating both of the frequency of the direct wave 
and the center frequency of the Doppler spread, the last frequency offset is estimated. Therefore, a sufficient frequency 
offset compensating range can be secured, and excellent BER characteristics can be realized. That is, securing of a 
sufficient frequency offset compensating range and realization of excellent BER characteristics are compatible with 
each other. Therefore, even under a Rician fading environment, compensation for fading distortion can be satisfactorily 
carried out. Therefore, the original data can be restored with high quality. Accordingly, service for receiver users can 
be improved. 

[0069] Fig. 10 is a graph showing the frequency offset estimated characteristics relating to this Embodiment 1 . In 
this figure, the circles show the BER in the case where the value obtained by synthesizing the frequency of the direct 
wave and the center frequency of the Doppler spread is regarded as frequency offset as in Embodiment 1 . The triangles 
show the BER when only the frequency of the direct wave is regarded as frequency offset. Furthermore, the squares 
show the BER when only the center frequency of the Doppler spread is regarded as frequency offset. 
[0070] As clearly understood from this graph, the BER characteristics relating to Embodiment 1 are more excellent 
than those in the case where only the frequency of the direct wave is regarded as frequency offset, and a frequency 
offset compensating range, which is wider than that in the case where only the center frequency of the Doppler spread 
is regarded as frequency offset, is realized in Embodiment 1 . 

[0071] Furthermore, according to the abovementioned Embodiment 1 , frequency offset candidates nAf RES are set 
for each estimation accuracy Af RES , signal powers P f (n) corresponding to the respective frequency offset candidates 
nAf RES are determined, and the signal powers P f (n) are smoothed by using the frequency windows, whereby the second 
frequency offset is estimated. Therefore, for example, if the estimation accuracies Af RES are finely divided, the estima- 
tion accuracy for the second frequency offset can be improved. Accordingly, further excellent BER characteristics can 
be realized. 

Embodiment 2 

[0072] Fig. 11 is a flowchart for explanation of the first frequency offset estimation processing relating to Embodiment 
2 of the invention. 

[0073] In the abovementioned Embodiment 1 , the phase difference 6 EP1 (mLN F ) is determined as the first frequency 
offset from the average phase difference vector D EPA (mLN F ). On the other hand, in this Embodiment 2, the phase 
difference 6 EP1 (mLN F ) is determined as the first frequency offset by further averaging the average phase difference 
vectors D EPA (ml_N F ) between the receiving burst signals by using a forgetting factor A.. Thereby, the estimation accuracy 
for the frequency offset is improved. 

[0074] Concretely, the DSP 26 determines the phase difference vectors D EP (mLN F ) (step W1 ), averages the phase 
difference vectors D EP (mLN F ) to determine the average phase difference vector D EPA (mLN F ) (step W2), and then 
executes interburst average processing (step W3). in the interburst average processing, the average phase'difference 
vectors D EPA (mLN F ) are further averaged between the receiving burst signals. 

[0075] More concretely, the DSP 26 holds the average phase difference vector D EPA (mLN F ) at least until the next 
receiving burst signal is received. Based on the average phase difference vector D EPA (mLN F ) when the receiving burst 
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signal B is received and the average phase difference vector <D EPA (mLN F )> B . 1 when the one previous receiving burst 
signal (B-1 ) is received, the DSP 26 determines the average phase difference vector <D EPA (mLN F )> B by way of Formula 
(19) below. In Formula (1 9), X is a forgetting factor, and is a value of 0 or more and 1 or less (0<X^1 ). 

5 <D EPA (mLN F )) B = D EPA (mLN F ) + MD E p A (mLN F )) B .., (1 9) 

[0076] Thereafter, the DSP 26 executes phase difference operation processing based on the average phase differ- 
ence vector <D EPA (mLN F )> B (step W4). Concretely, as shown by Formula (20) below, the DSP 26 determines the 
10 phase difference 8 EP1 (mLN F ) based on the average phase difference vector <D EPA (mLN F )> B . Thus, the DSP 26 esti- 
mates the first frequency offset coressponding to the frequency of the direct wave. 

= , , M . f -1 lm[(D EPA (mLN F ))] 

9 EP1 (mLN F ) = tan Re[(D EPA (mLN F ))] (20) 



75 



[0077] As described above, according to this Embodiment 2, average processing is executed while the previous 
average phase difference vectors D EPA (mLN F ) are gradually forgotten by using the forgetting factor X. Therefore, even 
in the case where the fading condition changes with the times, average processing can be executed for the average 
20 phase difference, vectors D EPA (mLN F ) while following the time variation. Accordingly, a phase difference vector in con- 
formity with the condition of the transmission path can be obtained. Therefore, even under an environment where the 
C/N is low and the frequency offset changes with the times, the frequency offset can be estimated with high accuracy. 



25 



Embodiment 3 



[0078] Fig. 1 2 is a flowchart for explanation of the second frequency offset estimation processing relating to Embod- 
iment 3 of the invention. 

[0079] In the abovementioned Embodiment 1 , the maximum value E^n^*) of the signal power is detected among 
the window signal powers E f (n) determined by one window signal power operation processing, whereby the second 

30 frequency offset is estimated. On the other hand, in this Embodiment 3, after averaging the window signal powers E f 
(n) between the receiving burst signals by using the forgetting factor X, the maximum value is detected to estimate the 
second frequency offset. Thereby, the estimation accuracy of the frequency offset is improved. 
[0080] Concretely, the DSP 26 determines the average distortion amounts C EP (kN F ) (step X1 ), and applies interblock 
DFT processing to the average distortion amounts C EP (kN F ), and thereby determines the signal powers P f (n) (step 

55 x2), and furthermore, determines the window signal powers E f (n) from the signal powers P,(n)(step X3), and then 
executes window signal power average processing (step X4). 

[0081] More concretely, the DSP 26 holds the window signal power E f (n) at least until the next receiving burst signal 
is received. Based on the window signal power E f (n) when the receiving burst signal B is received and the average 
window signal power <E f (n)> B . 1 when the one previous receiving burst signal (B-1 ) is received, the DSP 26 determines 
40 the average window signal power <E f (n)> B by way of Formula (21 ) below. In Formula (21 ) below, X is a forgetting factor 
between 0 and 1 (0<X<1). 

<E f (n)> B = E f (n) + X<E f (n)> B . 1 (21) 

45 

[0082] As shown by Formula (22) below, the DSP 26 detects the maximum value <E f (n max )> B average window signal 
power among the determined average window signal powers <E f (n)> B (step X5). 



50 (E f <n MAX )) B « M a nfedi)),] ( 22 ) 

[0083] Thereafter, the DSP 26 estimates the phase difference 6 EP2 (mLN F ) corresponding to the possible frequency 
offset nAf RES corresponding to this maximum value <Ef(n MAX )> B as the second frequency offset (step X6). 
55 [0084] As described above, according to Embodiment 3, average processing is executed while the previous window 
signal powers are gradually forgotten by using the forgetting factor X. Therefore, even in the case where the fading 
condition changes with the times, average processing of the window signal powers can be executed while following 
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can iTobT^ """r Si9nal P ° WerS " ° 0My With the C ° nditl0nS ° f the transmission channel 

S^finl™ TJ f u 60 Und6r aP env,ronment wher * C/N is low and the frequency offset changes with the 
times, frequency offset can be estimated with high accuracy. 

5 Embodiment 4 

«» Kimi n tinn m ^ 0C,iment *' t i nfluences of noise are elimi " a ^ by averaging the window signal powers E,(n) whereby 
the est.mat on accuracy of the second frequency offset is improved. On the other hand, improvement in the esttS 
accuracy of the second frequency offset by eliminating influences of noise can also te achteve tTaveraaXe 
s.gna. powers P, n). Therefore, in this Embodiment 4, the signal powers P ( (n) are a^Tte^ntoT^l 
burst s.gna.s to el.minate influences of noise, whereby the estimation accuracy of the second S5Z£lSES 

'5 [0087] Concretely, the DSP 26 determines the average distortion amounts C EPA (kN F ) (step Y1 ) and aoDlies DFT 
processmg to the average distortion amounts C EPA ( k N F ) to determine signa, powe1sP f (n) step Y2) and men executl 
n!«Z7 TT pr ° CeSSin 9 < s,e P Y3 >- M °re concretely, the DSP 26 holds the signal power R^taSTSSS 
next rece.v,ng burs, s.gnal is received. Based on the signa. power P f (n) when the receiving burst signa B^s received 
and the average s.gnal power <P,(n)> B . 1 when the one previous receiving burst signal (B 1) is received the DSP 26 

~oT< uoZT nal power <p,(n)>B by way of Formu,a (23) biL - ,n FomJa &sxsi5Es 



< P f( n »B = P f (n) + MP,m B ., (23) 

L 00 fh^ JT**" 6 !; J" 6 °t P 26 determines the window signal powers E f (n) based on average signal powers <P,(n) 

l!l T Y4> ' th6n det6CtS the maximum va,ue E < window signal power (step YS^ and 

est.mates the phase drfference e EP2 (mLN F ) corresponding to the frequency offset candidate nJZVnZonLlt 
the maximum value E,(n max ) as the second frequency offset (step Y6) ciaaie n A t RES correspond.ng to 

Sp 0 ^P^tTZl!lT n9 » 0 T Emb ° diment 4 - ,he ave ' a °o Pressing is executed while the previous 
E °T gradua " v ^gotten by using the forgetting factor X. Therefore, even in the case where the 

he ZTJr T Cha A n9eS H Wlth , the timeS ' thG aVera9e " rocessi "9 ° f the «ton-l Powers can be executed w^le foMow g 
obtained tZ ^ S ' 9na ' C ° nf0rmrty Wi,h the conditions ° f the transmission channefcTbe 

^^^^^^^ iS - "* - — offset changes w«b the times, 

Embodiment 5 

5 0 o 0 f 9 the invent ' eXP ' anati ° n ° f f ° ffSet estima,i °" P rocessi "9 bating to Embodiment 

nef the phase' SSSTJSS D^nTh^ "J" diSt0rti ° n am ° UntS ° f the -ion chan- 

ts Embodiment if htrft- ° ep( * Nf) between the ad i acent k "own signals is determined. On the other hand in 
LnLn T ^ ? determinm 9 onl V P has o difference information that is the scalar quantity between the adiacent 

Sonata fs7eo z Tln tl nh!l h « ? * dlfferenCe in,ormation operation processing first between the known 

IS J^"^^^^ 1 ? . ,n, r ati0 " ° Perati0n P rocessi "9. P*»» difference information as the scalar 

<^ZtL S * ™T ^ 38 Sh ° Wn by F ° rmU,a < 24 > below - ,he DSP 26 inverts 
nh««fS? am ° U " tS c EPi( kN F«) of the transm.ss.on channel into phase information epprfkMp+i) and determines 

^on T™ '" f T a,, ° n A6Ep(kNF) betWee " adjacenl known sionals 'n formula (24) belowl shows uTelSS 
cat.on numbers wh.ch correspond one to one known signals, and is a value between 0 and ^2 (0< <^ 2) 
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N p -2 



Atf EP (kNp) = J (H EPi (kN F + i + 1) - 6 EPi (kN F + i) ) 

N p - 1 ito 

6 EPi (kN F + i) - tan 1 — E — 

Re[c EPi (kN F +i)J 



[0093] Thereafter, the DSP 26 executes average processing for averaging the phase difference information AG EP 
(kN F ) in one burst (stepZ2). Concretely, the DSP 26 collects phase difference information AG EP (kN F ) for the estimation 
period Tf rom the estimation time (m-1 )T to the estimation time mT. If the number of known signal blocks in the estimation 
period T is L, the DSP 26 obtains L of phase difference information AeEpKm-IJLNp+INp). Herein, 1 shows the indication 
numbers which correspond one to one known signal blocks, and is a value between 0 and (L-1) (0<1<L-1). 
[0094] Thereafter, as shown by Formula (25) below, the DSP 26 applies average processing to the collected L of 
phase difference information Ae EP ((m-1)LN F +IN F ), whereby the average phase difference e EP1 (mLN F ) is determined 
as the first frequency offset. 

8 BP1 (mLNp) = - 2 1a0ep( {m " 1)LNf + 1Np) (25) 

L 1-0 



[0095] As described above, according to this Embodiment 5, the average processing is executed by using only the 
phase difference information. Therefore, the processing can be simplified. 

Embodiment 6 

[0096] Fig. 1 5 is a flowchart for explanation of the first frequency offset estimation processing relating to Embodiment 
6 of the invention. 

[0097] In the abovementioned Embodiment 1, the distortion amounts c EPj (kN F -i-i) of the transmission channel are 
multiplied by the complex conjugate and summed, whereby the phase difference vector D EP (kN F ) between adjacent 
known signals is determined. On the other hand, in Embodiment 6, by applying DFT processing to the distortion amounts 
c Ep .(kN F -M) of the transmission channel, the phase difference vector D EP (kN F ) is determined. 

[0098] Concretely, the DSP 26 executes inter-symbol DFT processing (step Rt). More concretely, by applying the 
D FT processing to the distortion amounts CEp^kNpH-i) of the transmission channel, the DSP 26 determines signal electric 
powers P f1 (n) corresponding to the plurality of frequency offset candidates nAf RES . 

[0099] Still more concretely, as shown by Formula (26), within the frequency offset estimating range f DET1 between 
-R s /2 and R s /2, the DSP 26 rotates the phases of the distortion amounts c EPj ((m-1 )LN F +IN F +i) of the transmission 
channel by the phase amounts corresponding to the respective frequency offset candidates nAf RES . Thereafter, the 
DSP 26 vector-synthesizes the distortion amounts Cgp^m-IJLNp+INp +i) whose phases have been rotated, and de- 
termines signal powers P f1 (n) corresponding to the plurality of frequency offset candidates nAf RES within the frequency 
offset estimating range. In addition, in Formula (26) below, n is approximately a value between -R s /(2Af RES ) and R s / 
(2Af RES )- 



1-0 



n pT 2 / . 2rcn(lN F 

£ c EPi ( (m - 1JLN F + 1N F + i) exp - j 

i-0 1 



|_ 2*n(lN F + i)Af j 



(26) 



[0100] Then, the DSP 26 vector-synthesizes the distortion amounts CEpitm-IJLNp+INp +i) whose phases have been 
rotated, and obtains signal powers P f1 (n) corresponding to the respective plurality of frequency offset candidates nA- 

*res- 

[0101] Then, the DSP 26 executes window signal power operation processing (step R2). Concretely, as shown by 
Formula (27), the DSP 26 sums the signal powers P f1 (n) corresponding to the frequency offset candidates nAf RES in 
the window having a frequency width of 2W 1 Af RES> whereby the window signal powers E^n) are determined. 
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w 

En(n) - § p fi( n + k ) (27) 



[01 02] Thereafter, the DS P 26 executes maximum value detection processing to detect the maximum value E. In 1 

zs?™ indow si T al f r rs E,i(n) (step R3) - and executes p hase differe " ce °p-tion 

quent o'set * R4) ' lhe PhaS6 difference G ^< mLN F> is ^imated as the firsHre 



quency offset. 
[0103] 



As described above, according to Embodiment 6, frequency offset candidates nAf BCC are set for each esti 

TZ^Vl™ T a ' P ° WerS Pfl(n) corres P° ndin 9 to the respective frequency oftsefcandidlTnll at 
determined and the signal powers P f1 (n) are smoothed by using the frequency window, whereby the first frequency 

JS U^^^^ - — ™ are finely set, the estimation Z=J 



Embodiment 7 



of invention 6 * ' eXP ' anati ° n ° f ,he ° ffs * Nation processing relating to Embodiment 7 

[0105] In the a'bovementioned Embodiments 1 through 6, the distortion amounts ^(kN-i) of the transmission 
dffer TT^Ts ^e aSed ZT PhaSe differenCe 6EPl(mLNF) aS the fet '^encTlet, and ,h~ 

tmas^ns. quency offset based ° n the dis,ortion am ° un,s c - (kN - i > aft - bei "9 — 

thlSct w^!! 1 ; '^t f ,r !K UenCy 0ffS6t COrreS P° nds to the frequency f D of the direct wave, and the frequency f D of 
the direct wave ex.sts near the center frequency f M of the Doppler spread (see Fig. 9A). Therefore on supposition thai 
1ST! ra, : 9e f arOUnd ,requenc y f o of the direct wave is a frequency offset estimating range, ^e^KSES 
that the center frequency f M of the Doppler spread exists within the range expected 

^J^^S^^JSS offset estimatin9 ran9e for estimation of the second ~ 

f!e^eLc C u°eTSi h , e ^ 26 T? diS, ° rti ° n am ° UntS c ^< kN ^> of transmission channel (step Q1 ), and 
then executes the first frequency offset estimation processing and the second frequency offset estimation orocessino 

^2cTS ! he DSP 26 e r ates thefirst frequency <^bJ^SSZZ 

offset C EPi (kN F +i) of the transmission channel. The DSP 26 temporarily holds, for example, this first frequency 

[0109] in addition the DSP 26 executes the first frequency offset estimation processing and distortion amount av- 
amfunf r Z 9 T\ U , ttane ° USly <SteP ° 2) - Concre t e 'V. the DSP 26 applies average processing to the dSortfon 
STSK; the transmission channel to determine the average distortion amount cJZj. 
[0110] Thereafter, the DSP 26 executes DFT processing (steps Q3 and Q4). Concretely the DSP 26 aoolies DFT 

EE« i n) COrrespondm 9 10 lhe P^rallty of frequency offset candidates nAf RES respectivefy determ ' neS S ' 9na ' 
S) and TerZZT* T 2 ! S6tS 3 freqU6nCy ° ffSet GStimatin 9 «W bSUd on the first frequency offset (step 
Si JlH T P ° WerS ,<n) COrres P° ndin 9 to the plurality of possib.e frequency offsets nAf RES within the 
frequency offset estimating range respectively (step Q4) RES 

wBh m ^ RE t S r ° r eaCh est,mation ^curacy Af RES . Herein, the frequency offset candidate in accordance 

oanZJx:r e z t is show : by NAfREs - ,n addi,ion - ,he dsp 26 sets »• ™» «™« *»• ^^ss 

f-VuLJ o« e f er '^ S , h ° Wn f ^ FOmr,U ' a (28) be, ° W ' b8Sed ° n the P hase amounts corresponding to the respective 
frequency offset candidates Af RES within this set frequency offset estimation range, the DSP 26 rotates the phases ^ 
the^erage distortion amounts c Ep((m+ l)N F) (0s1sL . 1) . In Formula (28) be|ow 9 n ' is . ^.^T^^ 
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Pf(n) 



L-l 



2 c EP ( (in + l)N F )exp - j 



RES 



(28) 



l-o 



[0114] The DSP 26 vector-synthesizes the average distortion amounts c EP ((m+1 )N F ) whose phases have been ro- 
tated. Thereby, the DSP 26 determines signal powers P f (n) corresponding to the plurality of frequency offset candidates 
nAf RES respectively. 

10 [0115] Thereafter, the DSP 26 determines window signal power E f (n) by cyclically summing the signal powers P f (n) 
in a frequency window with a predetermined frequency width (step Q5). The DSP 26 determines the maximum value 
E f (n MAX ) of the window signal power (step Q6), and estimates the frequency offset candidate nAf RES of this maximum 
value E f (n MAX ) as the last frequency offset (step Q7). 

[0116] As described above, according to Embodiment 7, the frequency offset estimating range to estimate the fre- 
15 quency offset based on the first frequency offset is set. Therefore, the first frequency offset elimination processing and 
phase synthesis processing become unnecessary. Therefore, the processing can be simplified. 
[0117] In addition, since the first frequency offset is estimated in a relatively wide frequency offset estimating range, 
a sufficient frequency offset compensating range is secured. Also : since the center frequency of the Doppler spread 
is used to estimate the last frequency offset, excellent BER characteristics are secured. Therefore, the same effects 
20 as in Embodiment 1 such as securing of a sufficient frequency offset compensating range and prevention of deterio- 
ration in the BER characteristics can also be obtained. 



Embodiment 8 

25 [0118] Fig. 17 is a flowchart for explanation of demodulation processing relating to Embodiment 8 of the invention. 
[0119] In Embodiment 1, the frequency offset elimination processing in which the phase of the digital baseband 
signal is rotated is executed after filtering processing. On the other hand, in this Embodiment 8, the frequency offset 
elimination processing is executed previous to the filtering processing. 

[0120] Frequency offset of the receiving burst signal occurs due to the stability of the oscillation circuit of the trans- 
30 mitter 1 and receiver 10 as explained in the subsection "Discription of the Related Art" above. Particularly, in the case 
where low-cost oscillation circuits are used, since the frequency stability is low, great frequency offset occurs. In ad- 
dition, the frequency stability of the oscillation circuits is influenced by the surrounding environment such as temperature 
and changes in the power supply voltage. 

[0121] Therefore, there is every possibility that the frequency offset becomes greater in accordance with time elapse 
35 although the offset is slight at the beginning. In this case, the frequency bandwidth of the receiving burst signal may 
expand to a large frequency bandwidth in comparison with the cutoff frequency in the filtering processing. Therefore, 
a part of the receiving burst signal may be cut by the filtering processing. Accordingly, the data decision may not be 
sufficiently carried out. 

[0122] In view of the above circumstances, in this Embodiment 8, the frequency offset elimination processing is 
40 executed prior to the filtering processing. Concretely, in this Embodiment 8, the frequency offset elimination processing 
is executed prior to the filtering processing, and distortion amount detection processing, frequency offset estimation 
processing, and integration processing are executed simultaneously with sequential processing including fading dis- 
tortion compensation processing and data decision processing. 

[0123] More concretely, as shown in Fig. 17A, when digital baseband signals are supplied from the A/D converter 
45 circuits 25a and 25b, previous to the filtering processing, the DSP 26 executes frequency offset elimination processing 
(step P1). Still more concretely, based on the accumulated phase difference 8((m-1)LN F +1 Np+i) determined by the 
separately executed integration processing, the phases of the digital baseband signals are rotated in reverse. Thereby, 
digital baseband signals from which frequency offset has been eliminated are obtained. 

[0124] Thereafter, the DSP 26 applies filtering processing to the digital baseband signals (step P2) to eliminate noise 
so components and the like. Thereafter, the DSP 26 detects a digital baseband signal corresponding to the Nyquist point 
(step S3), eliminates fading distortion from the digital baseband signal (step P4), and then executes data decision 
processing (step P5). 

[0125] On the other hand, as shown in Fig. 17B, the DSP 26 decides whether or not the digital baseband signal 
corresponding to the Nyquist point has been detected (step N1 ). If the Nyquist point has been detected, based on the 
55 digital baseband signal corresponding to the Nyquist point, the DSP 26 detects the distortion amount c EPj (kN F +i) of 
the transmission channel (step N2). Thereafter, based on the distortion amount c EPi (kN F +i) of the transmission channel, 
the DSP 26 estimates the phase difference G s (mLN F ) corresponding to the frequency offset (step N3), and cyclically 
sums the phase differences 8 s (mLN F ), whereby an accumulated phase difference e((m-1)LN F +IN F +i) is determined 
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(step N4). The DSP 26 uses this determined accumulated phase difference 6((m-1 )LN F +INjH-i) for the frequency offset 
elimination processing in step P1. 

[01 26] As described above, according to this Embodiment B f the frequency offset elimination processing is executed 
before the filtering processing. Therefore, even if the frequency bandwidths of the receiving burst signals are greater 
5 than the cutoff frequency in the filtering processing, the frequency offset can be eliminated while a part of the receiving 
burst signals is not cut. Therefore, data decision can be satisfactorily carried out. 

Embodiment 9 

10 [0127] Fig. 18 is a flowchart for explanation of demodulation processing relating to Embodiment 9 of the invention. 
[01 28] In the abovementioned Embodiments 1 through 8, by directly rotating the phase of the digital baseband signal, 
the frequency offset Is compensated. On the other hand, in this Embodiment 9, the phase of the receiving burst signal 
is rotated by changing the frequency of the local oscillation signal generated at the oscillation circuit 22, whereby the 
frequency offset is compensated. 

*5 [0129] Concretely, in this Embodiment 9, digital signal processing by way of software is not applied to the digital 
baseband signal, but a voltage to be applied to the oscillation circuit 22 comprised of VCO is controlled to eliminate 
the frequency offset. 

[0130] That is, in this Embodiment 9, as shown in Fig. 18A, filtering processing (step M1), Nyquist point detection 
processing (step M2), fading distortion compensation processing (step M3), and data decision processing (step M4) 
so are executed, and simultaneously with these sequential processings, as shown in Fig. 1 8B, distortion amount detection 
processing, frequency offset estimation processing, and VCO control processing for controlling the voltage to be applied 
to the VCO 22 are executed. 

[0131] More concretely, the DSP 26 decides whether or not the digital baseband signal corresponding to the Nyquist 
point has been detected (step L1). In the case where the Nyquist point has been detected, the DSP 26 detects the 
25 distortion amount c E pj(kNF+i) of the transmission channel based on the digital baseband signal corresponding to the 
Nyquist point (step L2). Thereafter, the DSP 26 estimates frequency offset based on the distortion amount c EPj (kN F +i) 
of the transmission channel (step L3). Then, the DS P 26 supplies a voltage in accordance with the estimated frequency 
offset to the VCO 22 (step L4). 

[0132] Thereby, the VCO 22 oscillates a local oscillation signal having an oscillation frequency which has deviated 
30 by an amount according to the frequency offset from the original oscillation frequency. Therefore, the frequency con- 
verter circuit 21 outputs an analog baseband signal from which frequency offset has been eliminated. Thus, excellent 
frequency offset compensation can be realized. 

[0133] As described above, according to Embodiment 9, by adjusting the frequency of the local oscillation signal 
generated by the VCO 22, the frequency offset is eliminated from the receiving IF signal. That is, before the filtering 
35 processing by the DSP 26, the frequency offset is eliminated. Therefore, as in the abovementioned Embodiment 8, 
even if the frequency bandwidth of the receiving burst signal expands to be larger than the cutoff frequency in the 
filtering processing, frequency offset can be eliminated while a part of the receiving burst signal is not cut. Therefore, 
the data decision can be satisfactorily carried out. 

40 Embodiment 10 

[0134] Fig. 19 is a flowchart for explanation of demodulation processing relating to Embodiment 10 of the invention. 
[0135] In the abovementioned Embodiment 1 , the first frequency offset and the second frequency offset are synthe- 
sized to estimate the last frequency offset, and in accordance with this frequency offset, compensation for the frequency 
offset is carried out, whereby the frequency of the digital baseband signal is automatically controlled. On the other 
hand, in this Embodiment 10, frequency offset compensation in accordance with the first frequency offset is carried 
out for the digital baseband signals, and frequency offset compensation in accordance with the second frequency offset 
is carried out for digital baseband signals afterthe abovementioned compensation is carried out, whereby the frequency 
of the digital baseband signal is automatically controlled. 

[0136] Concretely, ihe DSP 26 applies predetermined filtering processing to the digital baseband signals (step K1), 
and then detects a digital baseband signal corresponding to Nyquist point (step K2). Next, the DSP 26 executes the 
first automatic frequency control processing (hereinafter, referred to as thef irst AFC processing) forthe digital baseband 
signal (step K3), and then executes the second automatic frequency control processing (hereinafter, referred to as the 
second AFC processing) for the digital baseband signal after being subjected to the first AFC processing (step K4). 
55 Then, the DSP 26 executes fading distortion compensation processing for the digital baseband signal after being 
subjected to the second AFC processing (step K5), and applies data decision processing to the signal (step K6), 
whereby data corresponding to the original information signal is restored. 

[0137] Fig. 20 is a flowchart for explanation of the first AFC processing. The DSP 26 detects the distortion amounts 
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c EP j(kNp+i) of the transmission channel for each symbol (step J1 ), and then executes the first frequency offset estima- 
tion processing (step J2). The first frequency offset estimation processing is the same as that in the abovementioned 
Embodiment 1 , and by this processing, the DSP 26 estimates the phase difference e EP1 (ml_N F ) as the first frequency 
offset corresponding to the phase rotation amount for one symbol. 

5 [0138] Thereafter, the DSP 26 executes integration processing for cyclically summing the phase differences 0 EP1 
(mLN F ) (step J3) to obtain an accumulated phase difference e^m-^LNp+INp+i). In this case, the DSP 26 obtains the 
phase difference 8-, ((m-1 )LN F +IN F +i) accumulated from the estimation time mT for each symbol period T s . 
[0139] Thereafter, the DSP 26 executes the first frequency offset elimination processing (step J4). Concretely, based 
on the accumulated phase difference e^m-IJLNp+INp+i), the DSP 26 rotates in reverse the phase of the digital base- 

10 band signal r(kN F +i) by an amount in accordance with the accumulated phase difference e^m-l )LN F +IN F +i). Thereby, 
the first frequency offset can be eliminated from the digital baseband signal r(kNp+i). Thus, the DSP 26 achieves 
automatic frequency control of the digital baseband signals r(kN F +i) based on the first frequency offset. 
[0140] Fig. 21 is a flowchart for explanation of the second AFC processing. Based on the known signal blocks of the 
digital baseband signal after being subjected to the first AFC processing, the DSP 26 detects the distortion amount 

15 c E pi 2 (kNp+i) of the transmission channel for each symbol in the known signal blocks (step 11). 

[0141] Thereafter, based on the detected distortion amount c EPi2 (kN F +i) of the transmission channel, the DSP 26 
executes the second frequency offset estimation processing that is the same as in Embodiment 1 (step 12). Thereby, 
the DSP 26 estimates the phase difference e EP2 (mLN F ) as the second frequency offset corresponding to the phase 
rotation amount for one symbol. 

20 [0142] Thereafter, the DSP 26 executes integration processing for cyclically summing the phase differences 6 EP2 
(mLN F ) (step 13), and obtains an accumulated phase difference e 2 ((m-1)LN F +IN F +i). In this case, the DSP 26 obtains 
the phase difference 6 2 ((m-1 )LNp+1Np+i) accumulated from the estimation time mT for each symbol period T s . 
[0143] Thereafter, the DSP 26 executes the second frequency offset elimination processing (step 14). Concretely, 
based on the accumulated phase difference G 2 ((m-1)LNp+INp+i), the DSP 26 rotates in reverse the phase of the digital 

25 baseband signal r R1 (kN F +i) by an amount according to the accumulated phase difference G 2 ((nv1 )LN F +IN F -i-i). Thereby, 
the second frequency offset can be eliminated from the digital baseband signal rp^kNp-H). Thus, the DSP 26 achieves 
automatic frequency control of the digital baseband signal r R1 (kNp+i) based on the second frequency offset. 
[0144] As described above, according to this Embodiment 1 0, compatibility of securing of a sufficient frequency offset 
compensating range and securing of excellent BER characteristics is also achieved. 

30 

Embodiment 11 

[0145] Fig. 22 is a flowchart for explanation of demodulation processing relating to Embodiment 11 of the invention. 
[0146] In the abovementioned Embodiment 10, the frequency offset elimination processing is executed after the 
filtering processing. On the other hand, in this Embodiment 11 , the frequency offset elimination processing is executed 
before the filtering processing. 

[0147] Concretely, in this Embodiment 11 , the first frequency offset elimination processing that is a part of the first 
AFC processing is executed before the filtering processing, and the remaining processing of the first AFC processing, 
that is, distortion amount detection processing, the first frequency offset estimation processing, and integration process- 
ing are executed simultaneously with sequential processing. 

[0148] More concretely, as shown in Fig. 22A, when the digital baseband signals are supplied from the A/D converter 
circuits 25a and 25b, and previous to the filtering processing, the DSP 26 executes the first frequency offset elimination 
processing (step H1 ). Still more concretely, the DSP 26 rotates in reverse the phases of the digital baseband signals 
by amounts according to the accumulated phase difference e((m-1)LNp+IN F +i) corresponding to the first frequency 
offset determined by the integration processing executed simultaneously with the sequential processing. Thereby, 
digital baseband signals from which first frequency offset has been eliminated are obtained. 

[0149] Thereafter, the DSP 26 applies filtering processing to the digital baseband signals (step H2) to eliminate noise 
components and the like, and then detects a digital baseband signal corresponding to the Nyquist point (step H3). 
Thereafter, the DSP 26 executes the second AFC processing based on the digital baseband signal (step H4). By 
executing this second AFC processing, the DSP 26 can obtain a digital baseband signal in which the last frequency 
offset has been eliminated. Next, the DSP 26 eliminates fading distortion from the digital baseband signal (step H5), 
and then executes data decision processing (step H6). 

[0150] On the other hand, as shown in Fig. 22B, the DSP 26 decides whether or not the digital baseband signal 
corresponding to the Nyquist point has been detected (step G1 ). In the case where the Nyquist point has been detected, 
based on the digital baseband signal corresponding to the Nyquist point, the DSP 26 detects the distortion amount 
c EPi (kN F +i) of the transmission path (step G2). Thereafter, based on the distortion amount c EPI (kN F +i) of the transmis- 
sion channel, the DSP 26 estimates the phase difference e s (ml_N F ) corresponding to the first frequency offset (step 
G3), and cyclically sums the phase differences 6 s (mLN F ), whereby an accumulated phase difference e((m-t) 
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LN F +IN F +j) is determined (step G4). The DSP 26 uses this determined accumulated phase difference 9((m-1) 
LN f +IN f -h) for the first frequency offset elimination processing in step H1 . 

[0151] As described above, according to this Embodiment 11 , the first frequency offset elimination processing is 
executed before the filtering processing. Therefore, even when the frequency bandwidth of the receiving IF signal 
5 becomes larger than the cutoff frequency in the filtering processing, the first frequency offset can be eliminated while 
a part of the receiving IF signals is not cut. Therefore, data decision can be satisfactorily carried out. 

Embodiment 12 

10 [01 52] Fig. 23 is a flowchart for explanation of demodulation processing relating to Embodiment 1 2 of the invention. 
[0153] In the abovementioned Embodiment 10, by directly rotating the phase of the digital baseband signal, the 
frequency offset is compensated. On the other hand, in this Embodiment 12, by adjusting the frequency of the local 
oscillation signal generated by the oscillation circuit 22 that is the VCO, the phase of the receiving burst signal is rotated, 
whereby the frequency offset is compensated. 

[0154] Concretely, in this Embodiment 12, digital signal processing by way of software is not applied to the digital 
baseband signal, but the voltage to be applied to the oscillation circuit 22 formed of the VCO is controlled, whereby 
the frequency offset is eliminated. 

[0155] That is, in this Embodiment 12, as shown in Fig. 23A, the filtering processing (step F1), IMyquist point detection 
processing (step F2), the second AFC processing (step F3), fading distortion compensation processing (step F4) : and 
data decision processing (step F5) are executed, and simultaneously with these sequential processings, as shown in 
Fig. 23B, distortion amount detection processing, the first frequency offset estimation processing, and VCO control 
processing for controlling the voltage to be applied to the VCO 22 are executed. 

[0156] Concretely, the DSP 26 decides whether or not the digital baseband signal corresponding to the Nyquist point 
has been detected (step E1). In the case where the Nyquist point has been detected, based on the digital baseband 
signal corresponding to the Nyquist point, the DSP 26 detects the distortion amount c EPi (kN F +i) of the transmission 
channel (step E2). Thereafter, the DSP 26 estimates the first frequency offset based on the distortion amount c EPj 
(kNp+i) of the transmission channel (step E3). Then, the DSP 26 applies a voltage according to the estimated first 
frequency offset to the VCO 22 (step E4). 

[0157] Thereby, the VCO 22 oscillates the local oscillation signal having an oscillation frequency deviated by an 
amount in accordance with the first frequency offset from the original oscillation frequency. Therefore, the frequency 
converter circuit 21 outputs analog baseband signals in which the first frequency offset has been eliminated. Thus, 
excellent frequency offset compensation can be realized. 

[0158] As described above, according to this Embodiment 12, by adjusting the frequency of the local oscillation signal 
generated by the VCO 22, frequency offset is eliminated from the receiving IF signal. That is, the first frequency offset 
35 is eliminated before the filtering processing. Therefore, as in the abovementioned Embodiment 11, even when the 
frequency bandwidth of the receiving IF signal expands to be larger than the cutoff frequency in the filtering processing, 
frequency offset can be eliminated while a part of the receiving IF signal is not cut. Therefore, data decision can be 
satisfactorily carried out. 

40 Other Embodiments 

[0159] The Embodiments of the invention are as described above, however, the invention is not limited to the above- 
mentioned Embodiments. For example, in the Embodiments, the case where the TDMA is applied as the communica- 
tions method is described. However, as the communications method, the FDMA (Frequency Division Multiple Access), 
CDMA (Code Division Multiple Access), or the like may be applied. In such a case, the radio receiving signals are not 
burst signals but are continuous signals as in the case where TDMA is applied. However, needless to say, the invention 
can be easily applied even to this case. 

[0160] In addition, in the abovementioned Embodiments, the case where the demodulation processing is realized 
by way of software by the DSP 26 is described. However, of course, the steps of the demodulation processing to be 
50 executed by the DSP 26 may be realized by hardware circuits. 

Claims 
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An automatic frequency control method for controlling the frequency of a radio receiving signal by compensating 
frequency offset of the radio receiving signal periodically including a plurality of adjacent known signals, wherein 
the frequency of the direct wave of the radio receiving signal and the center frequency of the Doppler spread of 
the radio receiving signal are estimated based on the distortion amounts of the known signals included in the radio 
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receiving signal, and the frequency offset of the radio receiving signal is compensated based on both of these 
frequencies. 

An automatic frequency control method according to claim 1 , wherein the radio receiving signal is a burst signal 
synchronized with a predetermined time slot in the TDMA. 

An automatic frequency control method for controlling the frequency of a radio receiving signal by eliminating 
frequency offset from the radio receiving signal periodically including a plurality of adjacent known signals, com- 
prising: 

a frequency offset estimating step for estimating frequency offset of the radio receiving signal from the fre- 
quency of the direct wave of the radio receiving signal, which is estimated based on the distortion amounts of 
the known signals included in the radio receiving signal, and the center frequency of the Doppler spread of 
the radio receiving signal; and 

a frequency offset eliminating step for eliminating the estimated frequency offset from the radio receiving signal. 

An automatic frequency control method according to claim 3, wherein the frequency offset estimating step includes: 

a distortion amount operating step tor-determining the distortion amounts of the known signals included in the 
radio receiving signal; 

a first frequency offset estimating step for estimating the frequency of the direct wave of the radio receiving 
signal based on the distortion amounts between adjacent known signals among these determined distortion 
amounts; 

a distortion amount frequency offset eliminating step for eliminating frequency offset from the distortion 
amounts of the known signals based on the estimated frequency of the direct wave; 

a second frequency offset estimating step for estimating the center frequency of the Doppler spread of the 
radio receiving signal based on the distortion amounts between the periodically inserted known signals among 
the distortion amounts from which the frequency offset has been eliminated; and 

a synthesizing step for synthesizing the phase amount corresponding to the frequency of the direct wave and 
the phase amount corresponding to the center frequency of the Doppler broadening to estimate frequency 
offset of the radio receiving signal. 

An automatic frequency control method according to claim 4, wherein the first frequency offset estimating step 
includes steps of: 

determining the phase difference vectors based on the distortion amounts between adjacent known signals 
among the determined distortion amounts; 

determining the average phase difference vector by averaging the determined phase difference vectors for a 
predetermined period; and 

estimating the frequency of the direct wave based on the determined average phase difference vector. 

An automatic frequency control method according to claim 4, wherein the second frequency offset estimating step 
includes steps of: 

determining the average distortion amount by averaging the distortion amounts from which frequency offset 
has been eliminated; 

determining signal powers corresponding to a plurality of frequency offset candidates set at each predeter- 
mined interval within a predetermined frequency offset estimating range based on the determined average 
distortion amount; 

determining the window signal power of all frequency offset candidates within the frequency offset estimating 
range by summing the signal powers of the frequency offset candidates in the frequency window with a pre- 
determined frequency width among the determined signal powers; and 

estimating the frequency offset candidate corresponding to the maximum value of the determined window 
signal power as the center frequency of the Doppler spread. 

An automatic frequency control method according to claim 4, wherein the first frequency offset estimating step 
includes steps of: 
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determining phase difference information based on the distortion amounts between the adjacent known signals 
among the distortion amounts; 

determining average phase difference information by averaging the determined phase difference information 
for a predetermined period; and 

estimating the frequency of the direct wave based on the determined average phase difference information. 

An automatic frequency control method according to claim 4, wherein the first frequency offset estimating step 
includes steps of: 

determining signal powers corresponding to a plurality of frequency offset candidates set at each predeter- 
mined interval in a predetermined frequency offset estimating range based on the distortion amounts of the 
adjacent known signals among the distortion amounts; 

determining window signal powers corresponding to all frequency offset candidates by summing signal powers 
corresponding to the frequency offset candidates in a frequency windows with predetermined frequency widths 
among the determined signal powers; and 

estimating the frequency offset candidate corresponding to the maximum value of the determined window 
signal powers as the frequency of the direct wave. 

An automatic frequency control method according to claim 3, wherein the frequency offset estimating step includes: 

a distortion amount operating step for determining distortion amounts of the known signals included in the 
radio receiving signal; 

a first frequency offset estimating step for estimating the frequency of the direct wave of the radio receiving 
signal based on the distortion amounts between the adjacent known signals among the determined distortion 
amounts; 

an average distortion amount operating step for determining the average distortion amount by averaging the 
determined distortion amounts; 

a signal power operating step for determining signal powers corresponding to a plurality of frequency offset 
candidates set at each predetermined interval in a frequency offset estimating range regulated by the estimated 
frequency of the direct wave based on the determined average distortion amount; 

a window signal power operating step for determining window signal powers of all frequency offset candidates 
in the frequency offset estimating range by summing the signal powers of the frequency offset candidates in 
the frequency window with a predetermined frequency width among the determined signal powers; and 
a second frequency offset estimating step for estimating the frequency offset candidate corresponding to the 
maximum value of the determined window singal powers as the frequency offset of the radio receiving signal. 

An automatic frequency control method according to claim 3, further comprising a filtering step for eliminating high 
frequency components over the cutoff frequency from the radio receiving signal, wherein 

the frequency offset estimating step is for determining the distortion amounts of the known signals by using 
the radio receiving signal from which the high frequency components have been eliminated; and 
the frequency offset eliminating step is for eliminating the estimated frequency offset from the radio receiving 
signal whose high frequency components have not been eliminated by the filtering step. 

An automatic frequency control method for controlling the frequency of a radio receiving signal by eliminating 
frequency offset from the radio receiving signal periodically including a plurality of adjacent known signals com- 
prising: 

a first distortion amount operating step for determining the distortion amounts between the known signals 
included in ihe radio receiving signal; 

a first frequency offset estimating step for estimating the frequency of the direct wave of the radio receiving 
signal based on the distortion amounts of the adjacent known signals among the determined distortion 
amounts; 

a first frequency offset eliminating step for eliminating frequency offset corresponding to the estimated fre- 
quency of the direct wave from the radio receiving signal; 

a second distortion amount operating step for determining the distortion amounts of the known signals included 
in the radio receiving signal from which the frequency offset has been eliminated; 

a second frequency offset estimating step for estimating the center frequency of the Doppler spread of the 
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radio receiving signal based on the distortion amounts between the periodically included known signal blocks 
among the determined distortion amounts; 

a second frequency offset eliminating step for eliminating frequency offset corresponding to the estimated 
center frequency of the Doppler spread from the radio receiving signal. 

12. An automatic frequency control divice for controlling the frequency of a digital baseband signal by eliminating 
frequency offset from the digital baseband signal upon receiving the digital baseband signal outputted from an A/ 
D converter as an input corresponding to a radio receiving signal which periodically includes a plurality of adjacent 
known signals, wherein the frequency of the direct wave of the radio receiving signal and the center frequency of 
the Doppler spread of the radio receiving signal are estimated based on the distortion amounts of the known signals 
included in the digital baseband signal, and frequency offset is eliminated from the digital baseband signal based 
on both frequencies. 

13. A demodulator, comprising: 

a frequency converter circuit for converting a radio receiving signal periodically including a plurality of adjacent 
known signals into an analog baseband signal; 

an A/D converter circuit for converting this analog baseband signal into a digital baseband signal; 
a digital signal processing device : which receives the digital baseband signal generated by the A/D converter 
circuit as an input, estimates the frequency of the direct wave of the radio receiving signal and the center 
frequency of the Doppler spread of the radio receiving signal based on the distortion amounts of the known 
signals included in the inputted digital baseband signal, eliminates frequency offset from the digital baseband 
signal based on both frequencies, and eliminates fading distortion from the digital baseband signal from which 
the frequency offset has been eliminated and then demodulates the digital baseband signal. 

14. A demodulator according to claim 13, wherein 

the frequency converter circuit has a voltage control oscillation part, which oscillates a local oscillation signal 
forconverting the radio receiving signal into an analog baseband signal, and changes the oscillation frequency 
of the oscillation signal in accordance with an applied voltage; and 

the digital signal processing device eliminates frequency offset from the digital baseband signal by applying 
a voltage in accordance with both of the estimated frequencies to the voltage control oscillation part. 
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